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Abstract 

For more than 150 years the Navier-Stokes equations for thermodynamically quasi-equilibrium 
flows have been the cornerstone of modern computational fluid dynamics that underpins new fluid 
technologies. However, the applicable regime of the Navier-Stokes model in terms of the level of 
thermodynamic non-equilibrium in the local flowfield is not clear especially for hypersonic and 
low-speed micro/nano flows. Here, we re- visit the Navier-Stokes model in the framework of Boltz- 
mann statistics, and propose a new and more appropriate way of assessing non-equilibrium in the 
local flowfield, and the corresponding appropriateness of the Navier-Stokes model. Our theoretical 
analysis and numerical simulations confirm our proposed method. Through molecular dynamics 
simulations we reveal that the commonly-used Knudsen number, or a parametric combination of 
Knudsen and Mach numbers, may not be sufficient to accurately assess the departure of flowfields 
from equilibrium, and the applicability of the Navier-Stokes model. 
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The Navier-Stokes equations are the most widely-used model for fluid dynamics. Their 
impact is far-ranging, e.g. weather forecasting, modern transport system design including 
aeroplanes, cars, and ships, and energy generation from wind turbines. The fundamental 
assumptions for the Navier-Stokes model are that the fluid is a continuous medium, and the 
flow is close to thermodynamic equilibrium. 

For a gas, the Navier-Stokes equations can be regarded as a first order approximate 
solution to the Boltzmann equation, in terms of the Knudsen number Kn (the ratio of 
the molecular mean free path A to the system characteristic length L) [I]. The Boltzmann 
equation assumes both binary collisions between gas molecules and molecular chaos, so 
a single particle velocity distribution function can be employed. The Chapman-Enskog 
approach to the Boltzmann equation assumes 
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where the distribution functions can be asymptotically obtained from the Boltzmann 
equation. The Maxwell-Boltzmann equilibrium distribution, 

2RT 

is the zeroth-order solution f(°>, and leads to the Euler fluid equations. Here, p denotes the 
gas density, T the temperature, R the gas constant, and £ the peculiar velocity of molecules, 
which is £ — u where £ represents the molecular velocity and u is the macroscopic fluid 
velocity. In Eq.Q, provides a non-equilibrium correction of the order of Kn. To 
recover the Navier-Stokes equations, 
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where, p is the gas pressure, and the shear stress and the heat flux g« are related to the 
following first-order gradients of velocity and temperature, 



du<i dT 
a - = ~ 2 ^ q * = ~ K ^ (4) 



where /i and k denote the viscosity and thermal conductivity. We can continue the series to 
obtain a-order corrections to the distribution function in terms of the Knudsen number. 

The Chapman-Enskog technique indicates that the Euler equations are appropriate for 
thermodynamically equilibrium flows with Kn = 0, while the Navier-Stokes equations are 



valid for linear departures from equilibrium where the Knudsen number is close to zero. 
In high altitude applications, such as spacecraft re-entry into planetary atmospheres, or 
vacuum applications, low-pressure chemical processes, and micro/nano devices, the flows 
can be highly non-equilibrium and the Navier-Stokes equations fail to provide an adequate 
description. However, the Knudsen number alone may not be sufficient to describe the level 
of non-equilibrium in the flowfield, and to assess whether the Navier-Stokes equations are 
applicable or not. Both simulation data [2] and theoretical analysis [3] indicate that the 
level of non-equilibrium is also strongly influenced by the Mach number. 

Two types of breakdown criteria have often been used for indicating when the Navier- 
Stokes equations are appropriate ([21 SHZj and references therein). The first type is the 
local Knudsen number, X/(f)\d(j)/dx\, where <fi denotes a macroscopic flow quantity (typically 
density, temperature or pressure). The other type is the product of local Mach and Knudsen 
numbers, or its equivalent [2] . The first criterion indicates that the Knudsen number has to 
be small for the Navier-Stokes model to be valid. The second criterion requires the product 
of Ma and Kn to be small. However, the latter is not appropriate for the low Mach number 
flows typically occurring in micro/nano-devices, as the evidence is that the Navier-Stokes 
equations fail for a small KnMa but moderate Kn [8]. 

The inconsistency in these two types of criteria raises a fundamental question: is there a 
better way to assess the level of thermodynamic non-equilibrium in the local flowfield and 
the appropriateness of the Navier-Stokes equations? As the Navier-Stokes equations are the 
cornerstone of modern computational fluid dynamics (CFD), their validity range as a model 
should be clearly defined. Here, we address this fundamental issue, aiming to accurately 
assess the level of non-equilibrium of the local flowfield, and so redefine the applicability 
regime of the Navier-Stokes equations. 

Instead of using Knudsen and Mach numbers defined by macroscopic flow properties, 
the fundamental direct evidence for the level of thermodynamic non-equilibrium in the 
local flowfield resides in the distribution function itself. The distribution function may be 
considered in two parts: the equilibrium and the non-equilibrium components, i.e. 

f = r + r q . (5) 

In order to consider the validity of the Navier-Stokes equations, the non-equilibrium part 
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jneq ma y j^g f ur ther divided into two components, i.e. 

r* = + f {H) , (6) 

where is the first order non-equilibrium correction (at the Navier-Stokes level) which is 
given by Eq.([3|, and is the higher-order non-equilibrium correction beyond the Navier- 
Stokes level. So, the distribution function can be split into three components: 

f = f e q + fW + f (H) (7) 

If the non-equilibrium part f neq is negligible in comparison to f eg , then the flow is in equilib- 
rium, and the Euler equations can be recovered from the Boltzmann equation. Only when 
f( H ^ are the Navier-Stokes equations recovered from the Boltzmann equation. Using 
the distribution function directly we can thereby assess when the Navier-Stokes equations 
are valid, which is not only physically sound but also practical as many computational meth- 
ods provide information on the distribution function during simulations, e.g. direct solution 
of the Boltzmann equation [9], the lattice Boltzmann method [10J, the direct simulation 
Monte Carlo method [llj . and molecular dynamics [T3] . 

To evaluate how far the flowfield is away from equilibrium, we introduce a parameter Co 
to describe the departure from local equilibrium: 



Co 
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which is a relative error of / to the Maxwellian f eq . Similarly, a parameter C\ can be 
introduced to describe how far the flowfield is away from the Navier-Stokes regime: 



1 V/(/ (1) ) 2 ^ V /(/ (1) ) 2 ^ ' 

Here, C\ is a direct indicator of the relative error introduced by using a Navier-Stokes model 
on the flowfield. Together, Co and C\ provide both an accurate assessment of the level of 
non-equilibrium in the local flowfield, and an indication of the appropriateness of using the 
Navier-Stokes equations. 
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FIG. 1: Half-channel profiles of Co and C\ at Kn = 0.01, for various Ma. The top (Y = 0.5) and 
bottom (Y = —0.5) plates are moving with speeds U w in opposite directions. 



In the following section, we show quantitatively why the commonly used Knudsen and 
Mach numbers fail as local flowfield indicators, and how they are related to Co and C\. We 
use nonlinear shear-driven Couette flows as examples, where the two plates are moving with 
a speed of U w in opposite directions with their temperatures set to T . To obtain accurate gas 
molecular velocity distribution functions, we perform molecular dynamics (MD) simulations 
using the OpenFOAM code that includes the MD routines implemented by Macpherson 
and Reese et al. [12]. Monatomic Lennard- Jones argon molecules are simulated [13], and 
initially the molecules are spatially distributed in the domain of interest with a random 
Gaussian velocity distribution corresponding to an initially prescribed gas temperature. 
They are then allowed to relax through collisions until reaching a steady state before we 
take measurements. To achieve a smooth velocity distribution function, molecular velocity 
samples are then taken in every time step (0.001 r, where r = y^md 2 /e, with m being molar 
mass, d, the diameter of gas molecules, and e being related to the interaction strength of the 
molecules) for a total run time of at least 30000 r (in the extreme rarefied and high speed 
flow case below, up to 100000 r). We use 83500 molecules in each simulation, and assume 
diffuse gas molecule/wall interactions. 

The profiles of Co and C\ for various Couette flow cases are presented in Figures [T]j4| It 
is clearly shown that the level of non-equilibrium depends on both the Knudsen number and 
the Mach number {Ma = U w /\/RT ). Co, the indicator of departure from equilibrium, is 
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FIG. 2: Half-channel profiles of Co and C\ at Kn = 0.1, for various Mo. The top (Y = 0.5) and 
bottom (y = —0.5) plates are moving with speeds J7 TO in opposite directions. 
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FIG. 3: Half-channel profiles of Co and C\ at Kn = 0.5, for various Ma. 



small at Kn = 0.01 especially in the bulk region (Y=0). However, when the Mach number 
increases, Cq becomes larger. At a higher Mach number (Ma = 5), the values of Co are even 
similar to those in a flow with a much higher Knudsen number (see Kn = 0.5 and Ma = 0.2 
in Figure [3]), showing that the extent of departure from equilibrium can be similar for two 
flows with very different Knudsen numbers. The Knudsen number alone does not determine 
the level of thermodynamic non-equilibrium. 

For a Knudsen number of 0.01, the Navier-Stokes equations are usually regarded as valid. 
However, the value of Ci, which inversely indicates the appropriateness of using the Navier- 
Stokes equations, can increase with the Mach number, see Figure [TJ When the Mach number 
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FIG. 4: Half-channel profiles of Co and C\ at Kn = 1, for various Ma. The top (Y = 0.5) and 
bottom (Y = —0.5) plates are moving with speeds U w in opposite directions. 




is 1.0 and above, a significant proportion of the non-equilibrium flow information cannot 
be captured by the Navier-Stokes equations, even in this simple flow configuration. When 
Ma = 0.2, and in the bulk region (Y = 0), the Navier-Stokes equations hold because C\ <C 1. 
However, in the wall region, the Navier-Stokes equations are less accurate. The reason for 
this is the presence of Knudsen layer in the near-wall region, where the linear constitutive 
relation for stress that is assumed in the Navier-Stokes equations becomes inappropriate. 
When Ma increases to 1.0, higher-order fluid models may be required, even in the bulk, to 
accurately capture non-equilibrium information at conventional Knudsen numbers as low as 
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FIG. 6: Dependence of C\ on Kn and Ma in the bulk (left) and at the wall (right). 



0.01. 

For flows with Kn around 0.1, it is generally assumed that the Navier-Stokes equations 
will still be useful in the bulk flow region. Indeed, Figure [2] shows that less than 10% 
error will be introduced to the non-equilibrium distribution function if the Navier-Stokes 
equations are used for a flow with Ma = 0.2. However, C\ is larger than 0.1 in the bulk 
region when Ma is above 0.2, so substantial non-equilibrium flow information would not be 
captured by a Navier-Stokes analysis in these cases. Figure[3]also shows that non-equilibrium 
information cannot be properly captured, even for Mach numbers as low as 0.2, when the 
Knudsen number becomes large (e.g. 0.5). Figures [T|4] together show that neither the 
Knudsen number nor the simple product of Mach and Knudsen numbers can appropriately 
assess the level of thermodynamic non-equilibrium in flowfields. However, we can discover 
the appropriate dependencies of Co on Kn and Ma in certain flows. 

When the Knudsen number is small, f neq ^ as a first order approximation, hence, 
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For Couette flows, this can be simplified to 
Co 
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In the Navier-Stokes model, the velocity gradient turns out as 2U W /L over the whole flow- 
field. The temperature gradient can vary with position, but is zero at the centerline and 
~ 2U^h/kL at the wall. So we estimate C at the centerline from Eq.(ll) to be 



Co 



Lp pL 



Ma - KnMa } 



(12) 



and at the wall as 



Co 
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(13) 
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Therefore, the level of non-equilibrium increases with both the Knudsen and the Mach 
numbers. Note that this estimation is only appropriate when the Knudsen number is not 
too large. For large Kn, we need to analyse the relation numerically. For Ci, it is difficult 
to provide a solution even for small Knudsen numbers, as no generally agreed Burnett- 
order solution for is available. The numerical data shown in Figures [H] and [6] suggest a 
complicated general dependency of Co and C\ on the Knudsen and Mach numbers. 

In conclusion, we have addressed how to accurately assess the departure of flowfields 
from thermodynamic equilibrium, and how to identify when the Navier-Stokes model is 
applicable or not. Two new kinetic parameters based on the molecular velocity distribution 
function are proposed, with Cq assessing how far flowfields are away from equilibrium, while 
Ci indicating the validity of the Navier-Stokes model. Our MD numerical experiments 
conducted for Couette flows confirm that the flowfield and the appropriateness of the Navier- 
Stokes equations are not properly assessed by the Knudsen number alone, or even by the 
combination of Knudsen and Mach numbers proposed in the literature. The direct indicators 
Co and C\ show a complicated dependency on both Kn and Ma. Our proposed parameters 
are not only of theoretical interest, but could be used for model switching criteria in kinetic- 
continuum hybrid numerical schemes. 
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